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ABSTRACT 

Accurate knowledge of the non-linear dark-matter power spectrum is important for 
understanding the large-scale structure of the Universe, the statistics of dark-matter 
haloes and their evolution, and cosmological gravitational lensing. We analytically 
model the dark-matter power spectrum and its cross-power spectrum with dark-matter 
haloes. Our model extends the halo-model formalism, including realistic substructure 
population within individual dark-matter haloes and the scatter of the concentra- 
tion parameter at fixed halo mass. We consider three prescriptions for the mass- 
concentration relation and two for the substructure distribution in dark-matter haloes. 
We show that this extension of the halo model mainly increases the predicted power 
on the small scales, and is crucial for proper modeling the cosmological weak-lensing 
signal due to low-mass haloes. Our extended formalism shows how the halo model 
approach can be improved in accuracy as one increases the number of ingredients that 
are calibrated from n-body simulations. 
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1 INTRODUCTION 

Accurate studies of large-scale cosmic structures and galaxy 
clustering became possible with the advent of large galaxy 
redshift surveys (Guzzo et al. 2000; Meneux et al. 2009). Dif- 
ferent analyses have been carried out depending on galaxy 
luminosity, color and morphology (Norberg et al. 2002; Ze- 
havi et al. 2005; Padmanabhan et al. 2007). According to 
the standard scenario of structure formation, galaxies with 
dissimilar features reside in different dark matter haloes and 
have experienced varied formation histories. 

Galaxies are believed to form and reside in dark-matter 
haloes that extend much beyond their observable radii. Some 
of them are located at halo centers while others orbit around 
it, constituting the satellite population. Dark-matter haloes 
form by gravitational instability from dark-matter density 
fluctuations (Bond et al. 1991; Lacey & Cole 1993) and sub- 
sequently merge to form increasingly large haloes as cosmic 
time proceeds. Gas follows the dark-matter density pertur- 
bations. Once it reaches sufficiently high densities, dissipa- 
tive processes, shocks and cooling allow stars to form from 
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this gas (White & Rees 1978; Kauffmann et al. 1999). While 
the main lines of this scenario are widely accepted, its details 
are still poorly understood. 

The clustering strength of a given galaxy population 
is related to that of the dark matter halos which host the 
galaxies. It is possible to provide an accurate description 
of clustering in the small-scale non-linear regime even if one 
has no knowledge of how the haloes themselves are clustered 
(Sheth & Jain 1997; Smith et al. 2003). The halo-model of 
matter clustering (Scherrer & Bertschinger 1991; Peacock & 
Smith 2000; Seljak 2000; Scoccimarro et al. 2001; Cooray 
& Sheth 2002), which has been the subject of much recent 
interest, allows a parametrization of clustering even on large 
scales. 

To date, almost all analytic work based on the halo- 
model approach assumes that haloes are spherically sym- 
metric and that the matter density distribution around each 
halo centre is smooth. However, numerical simulations of 
hierarchical clustering have shown that haloes are neither 
spherically symmetric (Jing & Suto 2002; Allgood et al. 
2006; Hayashi et al. 2007) nor smooth (Moore et al. 1998; 
Springel et al. 2001; Gao et al. 2004; De Lucia et al. 2004; 
Tormen et al. 2004; Giocoli et al. 2008). About 10% of the 
mass in cluster-sized haloes is associated with subclumps. A 
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halo model which includes the effects of halo triaxiality on 
various clustering statistics is developed in Smith & Watts 
(2005); Smith et al. (2006), and formalism to include halo 
substructures was developed by Sheth & Jain (2003). 

The main purpose of this work is to include recent ad- 
vances in our understanding of halo substructures into the 
halo model approach. This is because an accurate model of 
substructures is a necessary first step to modeling the small 
scale weak gravitational lensing convergence and shear sig- 
nals (Bartelmann & Schneider 2001; Hagan et al. 2005) and 
for the substructure contribution to the gravitational flexion 
(Bacon et al. 2006). 

The present paper is organized as follows. In Sect. 2, we 
describe the ingredients of our extension of the halo model: 
the properties of subhaloes and the substructure mass func- 
tion. In Sect. 3, we show how to incorporate these into the 
halo model. Cross-correlations between haloes and mass as 
well as clumps and mass are studied in Sect. 4. Sect. 5 sum- 
marizes our methods and conclusions. 



2 THE MODEL 

We assume a flat ACDM cosmological model consistent with 
a combined analysis of the 2dFGRS (Colless et al. 2001) and 
first-year WMAP data (Spergel et al. 2003). It has the pa- 
rameters (fi m , Q b , h, n, ct 8 )=(0.25, 0.045, 0.73, 1, 0.9). This 
is the model adopted by Springel et al. (2005) for the Millen- 
nium Simulation (hereafter MS), against which we test our 
model. Our linear power spectrum uses the transfer func- 
tion by Eisenstein & Hu (1998), which takes baryonic acous- 
tic oscillations into account, in agreement with CMBFAST 
(Seljak & Zaldarriaga 1996) to 1%. 



2.1 Halos 

2.1.1 Mass function 

For modeling the non-linear dark-matter power spectrum, 
we require the redshift evolution of the number density of 
collapsed haloes, n{M,z). According to the spherical col- 
lapse model, a density fluctuation collapses when and where 
the linearly evolved density field smoothed with a top-hat 
filter exceeds the critical value 5 c (z) at that redshift. The 
distribution of collapsed peak heights is Gaussian when ex- 
pressed in terms of the variable fi(M,z) — S c (z) / \J S(M), 
where S(M) is the variance of a halo of mass M, if the col- 
lapse is not influenced by the surrounding gravitational field 
(Press & Schechter 1974; Bond et al. 1991; Lacey & Cole 
1993). Following this idea, one can write the halo number 
density in the form 
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where v = fi 2 (f(/j,)d(i is Gaussian for a fixed threshold S c ), 
and p is the mean matter density of the Universe. 

However, the dark-matter halo mass functions mea- 
sured in A-body simulations are far from Gaussian. The 
Press & Schechter (1974) formula over-predicts the abun- 
dance of low-mass haloes and under-predicts that of high- 
mass systems (Sheth & Tormen 1999). This is because the 
collapse of the haloes is affected by the gravitational tidal 



field of the surrounding matter (Sheth et al. 2001). Differ- 
ent functional forms have been proposed to fit the simula- 
tion results (Sheth & Tormen 1999; Jenkins et al. 2001). 
We use the Sheth & Tormen (1999) mass function because 
it is well tested against numerical simulations and easily 
parametrised in terms of the variable v which we shall use 
throughout. 

2.1.2 Bias 

Since the halo model assumes that all matter in the universe 
is in collapsed systems, we need to model the ratio between 
the power spectra of dark-matter haloes and of the dark 
matter, i.e. the bias parameter, to compute the dark-matter 
power spectrum. To be consistent with our choice of the 
mass function, we also use the bias parameter b{v) derived 
by Sheth & Tormen (1999). When expressed in terms of the 
peak height v, the bias parameter is independent of redshift. 
We recall that both the mass function and the bias parame- 
ter by Sheth & Tormen (1999) reproduce the original equa- 
tions by Press & Schechter (1974) and Mo & White (1996) 
when specialized to isolated spherical collapse. At fixed red- 
shift, the bias is an increasing function of the halo mass. 
More massive haloes (forming later and having a lower con- 
centration) are more biased compared to less massive haloes 
(forming earlier and with a higher concentration). This has 
been confirmed by recent A-Body simulations. 

Different improvements of the bias factor have been 
proposed (Jing 1998; Governato et al. 1999; Sheth et al. 
2001; Mandelbaum et al. 2005). Recent analyses have shown 
interesting features and unexpected results in view of the 
standard excursion-set formalism (Sheth & Tormen 2004), 
also confirmed by Gao et al. (2005); Faltenbacher & White 
(2010). At masses exceeding M* (i.e. peak heights v > 1), 
haloes in the same mass bin but with lower concentration 
(and thus lower formation redshift) are more biased than 
those with high concentrations, with a scatter of « 20% 
around the mean. The opposite holds for masses less than 
M* (y < 1). There is currently no accurate analytic model 
which allows a simple parametrization of the 'assembly bias' 
effect, so we have decided to continue using a bias parame- 
ter that depends on halo mass alone. Our algorithm can be 
easily extended by a new model for the bias factor if needed. 



2.1.3 Density profiles 

Various different definitions for the virial overdensity A^ 
have been used in the literature. Some authors (Jenkins 
et al. 2001; Gao et al. 2004) use A2000 corresponding to 200 
times the critical density, or A200&, defined as 200 times the 
background density (Diemand et al. 2007). Others (Sheth & 
Tormen 1999; Sheth et al. 2001; Bullock et al. 2001; Jing & 
Suto 2002) choose A^ = A v i r according to the spherical col- 
lapse model (Peebles 1980; Kitayama & Suto 1996; Bryan & 
Norman 1998). Here, we adopt A^ = A vir and use the virial 
overdensities given by Eke et al. (1996), which are related 
to the virial mass M v i r by 



AA^ZT 



47T 



(2) 



Qm(z) 

Cosmological iV-body simulations (Navarro et al. 1996; 
Power et al. 2003; Neto et al. 2007; Gao et al. 2008) have 
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Figure 1. Mass-concentration relations of dark-matter haloes at 
four different redshifts. The solid, dash-dotted and dashed curves 
show the models by Zhao et al. (2009), Neto et al. (2007) and 
Scljak (2000) at redshifts z=0, 0.5, 1 and 2, respectively. In the 
top-left panel, we also show the mass-concentration relation for 
substructures estimated by Pieri et al. (2009) fitting the results 
of the Aquarius simulation for clumps located at the virial radius 
(d v i r = 1) and at 0.1-R„i r from the host-halo centre (d v i r = 0.1). 
Notice that this fit has been rescaled to the definition of the virial 
radius adopted in this work. 



shown that the dark-matter density in isolated haloes has a 
universal radial profile whose slope steepens from —1 near 
the centre to —3 towards the virial radius, 



p(r\M vir 



r/r s (1 + r/r 3 ) 2 



(3) 



where r s is the scale radius of the halo. Its concentration is 
denned as c vir — R vir /r s . The amplitude p s sets the dark- 
matter density at the scale radius, 



R/Ivir 



ln(l + cw) 



C-vir 
1 ~\~ Cvi 



(4) 



Note that cw, and thus also p 3 , depend on halo mass; we 
expand on this in the next subsection. 

We shall need in the following an expression for the nor- 
malized Fourier transform of this density profile, truncated 
at the virial radius, 



u(k\M) 



4nr sin kr 
M kr 



p(r\M)dr 



(5) 



For the NFW profile, this integral has a convenient analytic 
expression (see Scoccimarro et al. 2001). 



2.1.4 Mass-concentration relations 

Motivated by numerical simulations and by previous halo- 
model studies, we consider three prescriptions for the mass- 
concentration relation to study uncertainties and differences 



in the recovery of the non-linear dark-matter power spec- 
trum by means of the halo model (Seljak 2000; Cooray & 
Sheth 2002). These are: 



(i) CI (Neto et al. 2007) 

N 1 \ 
Cvir \Z ) 



Cl4 



M 



(l+z) ^' 

This relation was obtained by Neto et al. (2007) by fitting 
the complete sample of haloes at z — in the Millennium 
Simulation (Springel et al. 2005). The value of the concentra- 
tion of a W 14 M Q /h halo estimated by Neto et al. (2007) for 
an enclosed mean density of 200 times the critical density 
has been rescaled to the virial overdensity. To rescale the 
normalization of the relation, we use the relation between 
C200 and Cvir presented in Maccio et al. (2008) at z = 0. 
(ii) C2 (Seljak 2000): 

-0.2 



(7) 



s M _ _9_ M 
vlrKZ) {l + z) [M, 

This steeper model was used by Seljak (2000) in order to fol- 
low the Peacock & Dodds (1996) fit for the standard CDM 
model adopted. The mass M*o represents the typical col- 
lapsed mass at z = such that <5 C (0) = <t(A/„o). Here, cr(M) 
is the rms fluctuation in the linear density-fluctuation field 
when smoothed with a top-hat filter which contains mass 
M. In this case, we set the concentration of an M* halo to 
9 in order to reproduce the non-linear power-spectrum of 
the Millennium Simulation at z = 0. Models CI and C2 in- 
corporate the redshift evolution of the mass-concentration 
relation proposed by Bullock et al. (2001). 
(iii) C3 (Zhao et al. 2009): 

., . , 8.41 1/8 

t(z) 



c z (z) 



1 + 
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(8) 



This model was proposed by Zhao et al. (2009) studying 
the mass-accretion history of dark-matter haloes in differ- 
ent cosmological models. They developed a model for the 
mass accretion history and related the concentration of an 
isolated halo to the time or redshift z when it had assem- 
bled 4% of its final mass. This is qualitatively motivated 
by the expectation that the structure of a halo should de- 
pend on its formation history. The correlation of the halo 
concentration with the time to. 04 is purely empirical. For 
a detailed explanation of how to implement model C3, see 
Appendix A of Zhao et al. (2009). So far, the model lacks a 
detailed physical justification. However, it expresses the no- 
tion that dark-matter haloes form first as consequence of a 
violent relaxation process (White 1996) which sets the value 
of the concentration near 4. Subsequent accretion adds ma- 
terial primarily to the outer shells. A minimum value of ~ 4 
means that the mass-concentration relation is expected to 
flatten at high redshifts, where haloes have had only had 
time to undergo the violent relaxation process. 

This flattening of the mass-concentration relation at high 
redshift was first discussed by Zhao et al. (2003b). Based on 
a combination of Af-body simulations with different resolu- 
tions, they found that the mass dependence of halo concen- 
trations weakens with increasing redshift. Later, Gao et al. 
(2008) and Duffy et al. (2008) confirmed these results study- 
ing haloes at different redshifts in the Millennium Simula- 
tion and with a different set of numerical simulations, re- 
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spectively. They proposed different power-law fits at various 
redshifts, but these are less complete than model C3. 

In Fig. 1, we show the mass-concentration relations de- 
scribed above at four different redshifts, z — 0, 0.5, 1 and 2. 
While models CI and C2 preserve their slopes as the redshift 
increases, model C3 flattens towards higher redshift. Also, 
notice that in model C3, no haloes have a concentration < 4, 
the value set by the violent relaxation process. 



2.1.5 Concentration distributions 

We shall adopt two ways of assigning concentrations to 
haloes of fixed mass, a deterministic model Dl and a stochas- 
tic model D2. 

(i) Dl: In common applications of the halo model, any 
scatter in the concentration parameter at fixed halo mass 
is neglected. The conditional distribution p(c\M) is thus as- 
sumed to be a delta distribution, 



p(c\M) = S[c-c(M,z}] 



(9) 



where c(M, z) is given by any of the mass-concentration re- 
lations CI, C2 or C3, see Eqs. (6, 7, 8). 

(ii) D2: However, numerical simulations have shown that 
at a fixed halo mass, different assembly histories yield dif- 
ferent values of the concentration (Navarro et al. 1996; 
Jing 2000; Wechsler et al. 2002). For example, Zhao et al. 
(2003a, b) proposed that the concentration could be related 
to the phase of fast accretion of the main halo progenitor. In 
general, haloes formed at higher redshift are more concen- 
trated than haloes formed more recently. The distribution 
around the mean value c found in numerical simulations is 
well described by a log-normal form, 



p{c\M) 



: exp 
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with variance oi nc . Values found in numerical simulations 
are 0.1 < a lnc < 0.25 (Jing 2000; Sheth & Tormen 2004; 
Dolag et al. 2004; Neto et al. 2007). 

Given a concentration distribution p(c\M), we define 
the marginalised density profile in Fourier space 



u(k\M) = / p(c\M)u(k\c(M))dc . 



(11) 



where c is again given by any of the Eqs. (6), (7) or (8). 
Note that u(k\M) is different from u(k\c(M)). 

2.2 Subhaloes 

2.2.1 Abundance and mass function 

Increasing the mass resolution in numerical simulations, re- 
cent studies have found that cores of accreted satellite haloes 
survive in host haloes as orbiting substructures (Tormen 
et al. 1997; Moore et al. 1999; Ghigna et al. 2000; Springel 
et al. 2001). Different studies have been performed in recent 
years, modeling the substructure mass function and its de- 
pendence on the mass and the redshift of the host halo (Gao 
et al. 2004; De Lucia et al. 2004; Giocoli et al. 2008; Angulo 
et al. 2009; Giocoli et al. 2010). Most of them find that the 
number of substructures of mass m, at fixed mass fraction 



m/M, depends on the host-halo mass: more massive haloes 
contain more subhaloes than less massive haloes. This has 
been translated by Gao et al. (2004) into the statement that 
the number of substructures per host-halo mass is univer- 
sal, recently also confirmed by Giocoli et al. (2010) using 
different post-processing of the same simulation. More mas- 
sive haloes assemble their mass quite recently, thus the time 
spent by their substructures under the influence of the grav- 
itational field of the host halo is shorter than in less massive 
hosts. Giocoli et al. (2010) showed that this also holds for 
haloes of similar mass, but different formation times. Sys- 
tems assembled more recently (thus with a lower concentra- 
tion) contain more substructures. 

In the model of Giocoli et al. (2010), the number of 
substructures per halo mass can be written as 



1 dN(c(M),z) 
M dlnm 



6N M (c(M) 
dlnm 
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(12) 



= (1 + z 



AMm a exp 



m) 



where Am = 9.33 x 10 -4 is a normalization factor and 
a = —0.9 and f3 — 12.2715 are the power-law slope and 
the steepness of the exponential cut-off. The term (1 + z) 1 ^ 2 
takes the redshift evolution of the normalization into ac- 
count. It depends on c(M tz )/c(M t o), while c/c reflects the 
dependence of the normalization on the concentration at 
fixed host-halo mass (Giocoli et al. 2010). Hereafter, subha- 
los within the virial region of their host will be called clumps. 
In the left panel of Fig. 2, we show the substructure mass 
function at z — in different host haloes. In the central 
panel we show the substructure mass fraction in terms of 
the host halo mass, 



m — dm , 

dm 



(13) 



where three different values of the minimum mass m m i„ are 
considered. In the right panel, we show the halo-occupation 
statistic of our model at z — 0, normalized by the number 
of subhaloes in a halo of mass M„, 

(N s ) = fV^dm, (14) 
Jo dm 

which is oc M and has {N S {N S - l))/(iV s ) 2 = 1. 
2.2.2 Spatial distribution 

Next, we need a model for the spatial distribution of sub- 
clumps around the host-halo centre. We study two spatial 
density distributions: 

(i) SI: The first assumes that the clumps follow the NFW 
dark matter distribution; 

(ii) S2: the second is motivated by numerical simulations, 
assuming that the subclump distribution is less concentrated 
than the dark matter (van den Bosch et al. 2004; Gao et al. 
2004). 

We shall later present results using the two models SI 
and S2 for the substructure distribution in the host to illus- 
trate by how much the power spectrum changes accordingly. 
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Figure 2. Left: Substructure mass function at the present time for different host-halo masses. Middle: Dependence of substructure 
mass fraction f s on host-halo mass. Curves show f a in subhalos above the minimum mass, as labeled. Right: Dependence of number of 
substructures on host halo mass, expressed in units of the number in an M*-halo. The solid line has a slope of unity. 



The normalized Fourier transform of the subclump distribu- 
tion in a spherically symmetric system is 

U s (k\c(M)) = f*™ ^l S -^L ns (r,c(M))dr. (15) 

For model SI, U s (k\M) = u{k\M). For model S2, we nu- 
merically transform the substructure number-density pro- 
file n s (r,M(c)) of Gao et al. (2004). In Fig. 3, we show the 
Fourier transform of the normalized density distribution for 
these two profiles. This shows that the Gao et al. (2004) 
fit drops at much smaller k than the NFW profile, because 
of the different trends of the two profiles towards the host- 
halo centre. We shall show in the next sections that this 
affects the non-linear power-spectrum at large k. The ring- 
ing in the Fourier transform is due to the fact that NFW is 
steeper than Gao et al. (2004) both in the center and at the 
virial radius. 



2.2.3 Mass-concentration relation 

The final ingredient we need to model is the mass- 
concentration relation for the substructures. In the top-left 
panel of Fig. 1, we show the mass-concentration relation fit 
to the substructures in the Aquarius simulation (Springel 
et al. 2008) by Pieri et al. (2009), 

c sub (m, dvir) = d~°: R [ci m^ 1 + c 2 m^ 2 ] ■ (16) 

Here, d V i r = r/R V i r is the distance from the host-halo 
centre in units of the virial radius. The parameters are 
listed in Tab. 2 of Pieri et al. (2009) for haloes 200 times 
denser than the background. In the figure, we plot the mass- 
concentration relation for clumps at 10% and 100% of the 
virial radius. This proposed relation takes into account that 
substructures with the same mass are more concentrated if 
they are closer to the host-halo centre, which is necessary 
for their survival. To be consistent with our choice of the 
virial overdensity, we have normalized this distribution to 
A vir . Since for subclumps at the virial radius (d vir — 1), 
the fit (16) follows the mass-concentration relation of iso- 
lated haloes and increases with decreasing distance from the 
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Figure 3. Fourier transform of the density profile. The dotted 
curves show the analytic Fourier transform of the NFW profile of 
different host-halo masses: 10 15 , 10 13 , 10 11 and 10 9 M o /h. The 
solid curves show the numerical Fourier transform of the Gao et al. 
(2004) substructure density distribution. The mass-concentration 
relation by Zhao et al. (2009) was adopted here. 

host-halo centre, we adopt for each model 

Csub (^l; d> v i r *) = d, v ^ r C v ^ r , (1^) 

where I =N, S, or Z, corresponding to CI, C2 and C3 (see 
Eqs. 6, 7, 8 and 16). 

2.3 Summary 

Before we proceed, let us briefly summarize the main ingre- 
dients of our extension to the halo model. 

• The abundance and bias of host haloes are taken from 
Sheth & Tormen (1999). 
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• Host haloes and subhaloes are assumed to have NFW 
density profiles. 

• We study deterministic (Dl) and stochastic (D2) mod- 
els of halo concentration at fixed mass, and explore three 
different choices for the mean relation between mass and 
concentration (CI, C2, C3). 

• We study two choices for the density run of subhaloes 
around the center of the host: either it is the same NFW 
profile as the smooth component (SI), or it has the form 
suggested by Gao et al. (2004) (S2). 

• The mass function of subhalos is given by equation (12), 
found in and normalized to numerical simulations. 

• The concentration of subhaloes increases towards the 
centre of their host halo, in agreement with numerical sim- 
ulations (Pieri et al. 2009). As we describe below, our halo 
model includes this effect only rather crudely. 



3 NON-LINEAR POWER SPECTRUM 

In this Section, we describe the halo-model approach at the 
non-linear dark-matter power spectrum and study how it 
changes at small scales when different models for the mass- 
concentration relation are taken into account. We will also 
show how the non-linear power spectrum can be decomposed 
when a realistic model for the substructure mass function in 
host haloes and two models for their spatial distribution 
are taken into account. The formalism for the substructure 
contribution to the one-halo term of the non-linear dark- 
matter power spectrum has been developed by Sheth & Jain 
(2003). Here, we shall also consider the clump contributions 
to the two-halo term and show the relevant equations. 

3.1 Contributions from host haloes 

In the halo model approach (Scherrer & Bertschinger 1991; 
Seljak 2000; Scoccimarro et al. 2001; Cooray & Sheth 2002), 
the two-point dark-matter correlation function is 

£(£ - x) = £ih(£- x) + &h{£- x) , (18) 

where the first or Poisson term describes the contribution 
to the matter density from individual haloes, while the sec- 
ond term describes the contribution from halo correlations. 
While both terms require knowledge of the halo mass func- 
tion and their dark-matter density profiles, the second needs 
also the two-point correlation function of haloes of different 
mass £hh(r\Mi, M2). Since the two-point correlation func- 
tion on large scales is dominated by the two-halo term and 
has to follow the linear correlation function, a good and con- 
venient way to express the two-halo correlation function is 

&fc(r|Mi, M a ) « b(Mi)b(M 2 )tu n (r) . (19) 

In view of later convolutions, it is convenient to work in 
Fourier space. We then define the dark matter power spec- 
trum P(k, z) as 

P(k, z) = 4tt f £(r, z) Sln ( fcr ) r 2 dr . (20) 
J kr 

For further convenience, and as is common in the large scale 
structure community, we define the dimensionless quantity 

A 2 (M) = ^%^. (21) 



The dark matter power spectrum can be decomposed as 

P(k,z) = P 1H (k,z) + P 2H (k,z) (22) 

where 

Pi H (k,z) = J (^j n(M,z) 

x j p{c\M)u 2 (k\c(M))dcdM, (23) 

and 

P2 H {k,z) = P lin (k)[j jn{M,z)b(M,z) 

x J p(c\M) u(k\c, M) ^m'z) dcdM ] 2 (24) 

are the one- and two-halo contributions. Here, b(c, M, z) is 
the bias factor with respect to the dark-matter of halos of 
mass M at z having concentration c. The dependence of 
halo bias on c as well as M is a simple way of including the 
'assembly bias' effect in the model. In practice, we ignore this 
effect, meaning that we set b(c, M, z)/b(M, z) = 1 from now 
onwards. Note that the integrals extend over all collapsed 
masses. 1 

The most commonly used approximations in the litera- 
ture (Seljak 2000; Cooray & Sheth 2002) use the determinis- 
tic model Dl for the concentration distribution and consider 
steep models for the mass-concentration relation in order to 
not loose small scale power. 

In Fig. 4 we show the non-linear dark-matter power 
spectrum predicted using the halo model, as described 
above, divided by the power-spectrum fit by Smith et al. 
(2003) at four different redshifts (from top to bottom) and 
assuming three models for the mass-concentration relation 
(from left to right). The solid curve refers to the analyti- 
cal non-linear power spectrum with the deterministic con- 
centration model (Dl), while the dotted and dashed curves 
refer to the model which adopt the stochastic model (D2) 
with widths a\ nc = 0.25 and 0.35, respectively. The ringing 
of the curves for large k is due to the Fourier tranform of 
the density profile. The thick solid curve refers to the power 
spectrum measured in the Millennium Simulation (Springel 
et al. 2005; Boylan-Kolchin et al. 2009). The dotted-dashed 
line in the right panels shows the shot-noise limit in the 
Millennium Simulation, which we do not remove from the 
numerical power spectrum, at the corresponding redshifts. 

From the Figure, we notice that model CI reproduces 
the MS power spectrum quite well up to k w 10, while model 
C2 performs better. Model C3 is also quite accurate and does 
not decrease at small scales, although it somewhat overes- 
timates the power in the range 1 < k < 50. Note that the 
stochastic concentration model (D2) increases the power on 
small scales (k > 10) by 15 to 20% compared to the deter- 
ministic (Dl) case. 

3.2 Contributions from subhaloes 

Considering the contribution from substructures, the 1-halo 
term can be further decomposed into four contributions from 

A lower bound of zero in the mass integrals actually means that 
the integration starts at 10 — 6 Mq (Green et al. 2004; Giocoli et al. 
2008). 
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Figure 4. Ratio of the non-linear power spectrum predicted by the halo model and the fit by Smith ct al. (2003). From left to right, 
the panels show results for the three models of the mass-concentration relation, as described in the text. From top to bottom, we show 
results at four different redshifts. In all panels, solid curves show the ratio for the deterministic concentration assignment, while dotted 
and dashed curves refer to the stochastic model with widths <ri nc = 0.25 and 0.35. The heavy solid curves show the ratio between the 
power spectrum in the Millennium Simulation and Smith ct al. (2003) at the corresponding redshifts, and the dash-dotted line (only in 
the right panels) shows the shot noise in the MS with respect to the fit. 
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the mutual correlations between smooth and clump compo- 
nents as explained in what follows. 



(1) Smooth-smooth correlation: 

PlH,ss(k,z) = 



j n(M,z) 



(25) 



u 2 (k\c(M sm ))p(c\M)dcdM 

(2) Smooth-clump correlation: 

PiH, sa {k,z) = 2 J (Jj) n(M,z) f ^u(k\c(M sm )) 
x Us(k\c(M))X c(M) (k,z)p(c\M)dcdM (26) 

(3) Two-point correlation between different clumps in the 
same halo: 

PiB.cc(k,z) = f (j) n(M,z) 



(27) 



x j U 2 s {k\c(M))i 2 c{M) {k,z)p{c\M)dcdM 



(4) Two-point correlation between pairs in the same clump: 
PiH,sei f - c (k,z) = f (jfj n(M,z) (28) 
x j J c(M) (k,z)p(c\M)dcdM 



Here, we define the smooth mass as 
M sm = (1 - f s ) M, 



(29) 



which depends on the host halo concentration (see eq. 12), 
while X and J are the functions 

^ f M m f . dN(c(M), z) , 
l c(M )(k,z) = j —j u m (k,c) ; ' > p(c\M)dcdm 



dm 



(30) 



and 



J cW (k,z) = ^yfuUk,c) dN ^ Z \ (clM)dcdrr< 

(31) 

where we also take into account the scatter in at fixed 
mass in the subhalo mass-concentration relation when the 
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Figure 5. Top: Ratio of the 1-halo components divided into substructures (Eqs. 26, 26, 28 and 29) compared to the 1-halo term (Eq. 23). 
A width of <T[ n c = 0.25 was assumed in the concentration distribution. In the left and right panels, we adopt the NFW and the Gao et al. 
(2004) model for the substructure distribution. Bottom: Ratio of the 2-halo components (Eqs. 32, 33 and 34) compared to the 2-halo 
term (Eq. 24). 



stochastic model is considered. The subhalo concentration 
depends on its mass and also on its radial distance from the 
host halo centre (see eq.17). The mass-concentration rela- 
tion adopted is consistent with the one assumed for isolated 
host haloes, such that the subhalo mass-concentration rela- 
tion, for distance larger then the virial radius of the host 
follows the halo mass-concentration relation. To account for 
the radial distance dependance we randomly sample the sub- 
structure density profile distribution in the host, consistently 
with the model adopted (NFW or Gao et al. (2004)), and 
locate the subhalo at a distance d V i r from the host centre 
(d v i r is expressed in unit of the host halo virial radius). Note 
that we do not consider scatter in / s at fixed c and M. 

Likewise, the large-scale (2-Halo) term can be further 
decomposed according to the different correlations between 
smooth and clump components. In this case, we have the 
three contributions 

(1) Smooth-smooth component correlation on large scales: 



P2H,ss{k,z) = Pu n (k)S l2 (k,z) 



(32) 



(2) Smooth-clump correlation: 

P2H,sc(k, z) = 2P Un {k)S I {k, z)C x (k, z) (33) 



(3) Clump correlation: 



P 2 H,cc(k,z) = P Un (k)C l2 {k, Z ) 



Here, 



M 



(34) 



S J (fc, 2 ) = / —n(M,z)b(M,z) 



J ^ u ( fc > c ( M s"0)P( c l M ) dcdM > (35) 



and 



M 



C'faz) = / —n(M,z)b(M,z) 
P 



x J I c{ M)(k,z)U s (k\c(M))p(c\M)dcdM. (36) 

Recall that / = N,S, or Z, corresponding to CI, C2 and C3 
(see Eqs. 6, 7, 8 and 16). 
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Figure 6. Non-linear dark-matter power spectrum, with substructures, at four different redshifts. We show three models referring 
to different mass-concentration relations. The heavy solid line shows the data from the Millennium Simulation. The width of the 
concentration distribution is assumed to be tri nc = 0.25, and the substructure distribution in the host halo is modeled according to Gao 
et al. (2004). For all models, the mass distribution in both the host haloes and in their substructures follows the NFW profile. 



The full non-linear dark-matter power spectrum is the 
sum of all components, 

4 3 

P+(k,z)=^2PiH,i(k,z)+^2P2H,i(k,z). (37) 

i=i i=i 

In Fig. 5 we show the ratio of the individual contribu- 
tions to the 1- and 2-halo components. We set the width 
fine = 0.25 for the concentration distribution, and adopt 
C3 for the mass-concentration relations. All functions de- 
pending on the concentration are estimated by marginaliz- 
ing over c. Models SI and S2 for the subclump distribution 
within the host are adopted in the different panels. Note 
that in both the 1- and 2-halo terms, the dominant contri- 
bution is due to pairs in which both members are in the 
smooth component; smooth-clump and clump-clump pairs 
never contribute more than 10% of the signal, and they be- 
come even less dominant at large k. The transition scale is 
around k w 1, since it represents the size of a typical col- 
lapsed halo in which the clumps are located. For k > 1 the 
smooth-clump and clump-clump terms start to drop down 
due to the density profile distribution of subhaloes in the 
host system. In other words, for scale smaller then 1 Mpc 



ft -1 the two-point probability has a threshold due to the 
typical distribution scale of the clumps within the virial ra- 
dius of the host 2 . Comparing right and left panels we notice 
that the power-spectrum components are quite sensitive to 
the more rapid decrease of the Fourier transform model S2 
(right panels) compared to model SI (left panels), both in 
the 1-Halo and the 2-Halo terms. However, the contribution 
from pairs within the same subclump (self-clump) increases 
with k reaching values of order 50% at k w 10 4 , this is due to 
the matter density profile distribution in small clumps, more 
and more concentrated as their mass decrises. The probabil- 
ity to find a pair of points of the same clump becames larger 
and larger as the scale decreases - so for large value of k. 

In Fig. 6, we show the full non-linear power spectrum at 
four different redshifts. We include results for models CI, C2 
and C3 and for a\ nc = 0.25. We adopt the NFW profile for 
the matter-density distribution in both the host haloes and 



2 We recall that the power spectrum is the Fourier tranform of 
the two-point correlation function. 
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Figure 7. Ratio between the predicted power-spectrum in our ex- 
tension of the halo model, which includes substructures and the 
stochastic assignment of halo concentrations, to the power spec- 
trum in the standard implementation of the halo model, which 
ignores halo substructures and assigns concentrations determin- 
istically. For the log-normal scatter in concentration we adopt 
°"lnc = 0.25. The solid and dotted curves assume an NFW or 
a Gao et al. (2004) model for the subclump distribution. The 
dashed line shows the predicted ratio when we allow scatter in 
the concentrations but assume smooth halos with no substruc- 
ture. The thick magenta line shows the power-spectrum ratio for 
the Millennium Simulation and the dash-dotted line shows the 
shot-noise. 



their substructures, and model S2 for the radial subclump 
distribution within the host haloes. 

In Fig. 7 we show the ratio between the full non-linear 
power spectrum with <ri nc = 0.25, considering models SI 
(solid) and S2 (dotted), and the standard halo model - 1 and 
2-Halo terms and no scatter in concentration, at four differ- 
ent redshifts as in Figure 6. This figure quantifies the impor- 
tance of substructures and the concentration scatter in the 
halo model. Including both effects increases the small-scale 
power by « 30%. On the large scales where the 2-Halo term 
dominates, accounting for or ignoring substructure makes 
no difference, as expected (Sheth & Jain 2003). In the figure 
we also show the ratio between the power spectra in the MS 
(thick solid line) and the standard halo model with stochas- 
tic concentration assignment (dashed line), compared to 
the standard halo model without concentration scatter. The 
dash-dotted line show the shot-noise. 



We shall also quantify the scales on which different host- 
halo masses and substructures contribute significantly. In 
contrast to Hayashi & White (2008), who presented a sim- 
ple model for the auto- and cross-correlation functions, we 
study all contributions to the non-linear power spectrum ex- 
plicitly. As for the power-spectrum in the previous section, 
our model of the 1-halo term follows Sheth & Jain (2003). 



4.1 Cross-correlation of haloes and mass 

We begin by writing down the contributions to the 1-halo 
term of the cross-correlation power spectrum between haloes 
and matter. Sitting at the center of a host halo, and as- 
suming that the matter is divided into smooth and clumpy 
components, the halo-smooth self-correlation and the halo- 
clump cross-correlation can be written as 



^ 1H,sH>(^J Z ) ~ 
X 

and 



M 



/M 
-^u{k\c{M sm ))p(c\M)dcdM, (38) 



— n(M,z) 
P 



x / U a (k,c)l c{M) (k,z)p(c\M)dcdM. (39) 

We recall that M 3rn is the smooth mass of the host halo, 
u(k\c(M sm )) and U s (k\c(M)) are the Fourier transforms of 
the density profiles of the smooth component and of the 
clump distribution respectively. Moreover, n(M, z) is the 
halo mass function at redshift z and I c (m) is given by 
Eq. (30). 

Regarding the contributions to large scales (the 2-halo 
term), we start from the same idea, i.e. we imagine sitting at 
the center of a host halo and cross-correlate with the smooth 
and the clumpy components of a distant halo. However, in 
this case we also have to take the halo bias relative to the 
dark matter distribution into account, as well as the non- 
linear dark-matter power spectrum. The smooth and clump 
cross-correlation power spectra are 



i&.H.CM = P ""i M) / n(M,z)b(M,z)dMS 1 (k,z) 

(40) 



rih 



and 



= PuJjt, z) r n , MjZ } h , Mz ^ dMC i, ktZ \ 

fih J 



with S 1 and C 1 given by Eqs. (35) and (36). 



(41) 



4 CROSS CORRELATIONS 

We now describe halo models of the cross correlation be- 
tween haloes and mass, as well as between substructures 
and mass. We shall show that while the cross correlation 
with the smooth component dominates the signal on large 
scales, the self-cross correlation dominates on small scales. 



4.2 Cross-correlation of subhaloes and mass 

Let us now place ourselves at the center of a clump and de- 
termine the cross-correlation on small scales. In this case, 
the 1-halo term will be the sum of three components: (i) the 
self-correlation with the substructure mass; (ii) the cross- 
correlation with the mass in the other substructures con- 
tained in the same host halo, and (iii) the cross-correlation 
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Figure 8. Halo and subhalo-mass cross-correlations at redshift z = 0.1 (top) and z = 1 (bottom). In each panel, we show the contributions 
of haloes and substructures to the 1- and 2-halo terms. The left and right panels show the cross-correlation between haloes and mass and 
between substructures and mass, respectively. We adopt the Gao et al. (2004) model for the radial subclump distribution in the host 
halo, Zhao ct al. (2009) for the mass-concentration relation and &\ nc = 0.25 for its width. The heavy solid curve shows the sum of all 
contributions. For comparison, we also show in each panel the total cross-correlation function when the spatial distribution of subclumps 
in the host follows the NFW model (heavy-dotted curve). 



with the smooth component. These three terms are 

P? Bt8m (k,z) = i| jn(M,z) J ^u(k\c(M am )) 

x Us (k\c(M))N s (c(M) , z)p{c\M)dc dM(42) 

further 

P? BlSe {k,*) = - / -n(M,z) f U 2 s {k\c{M)) 
n s J p J 

x X c(M) (k, z)N s {c(M), z)p(c\M)dcdM (43) 



and 
P 



iSr,.Sc(*.*) = j —n(M,z) J X c(M) (k,z)p(c\M)dcdM, 

P (44) 
where we have underlined that N s , the number of substruc- 
tures in a M-halo, depends on the mass through the con- 
centration and on redshift. For the 2-halo term, we have 
two equations analogous to those of the halo-matter cross- 
correlation, and we need to integrate over the substructure 



mass function. We find 

Pun{k,z) 



,{k,z) = 



j n(M,z)b{M,z) J N 3 (c(M),z) 
x Us(k\c(M))p(c\M)dcdMS I {k,z), (45) 



and 



2H,S, 



(k,z) 



n{M,z)b{M,z) I N a (c(M),z) 



Plin{k, Z) 



x (7 s (fe|c(M))p(c|M)dcdMC 7 (fc,2). (46) 



We show the results in real space for each component of the 
cross-correlation after the inverse Fourier transform 



t i \ 1 f D n N sin(fcr) 
6,x(r,*) = / P lH , k {Kz)—^— 



(47) 



where the index i runs over all halo and substructure cross- 
correlation terms. Going from Fourier to real space, some of 
the 1-halo terms are simple, because we are inverse Fourier 
transforming the NFW density profile and the density run 
of subclumps in the host (see eq. 15). For the model SI, in 
which subclumps follow the dark matter distribution, terms 
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Figure 9. Halo and subhalo-mass cross-correlations at z = 0.1 and 2 = 1 separating the contributions from different host-halo masses. 
The left panels show the halo-mass cross-correlation. The 2-halo terms in different mass bins have been normalized to the total number of 
haloes between 10 10 and 10 14 Mq//i. The right panels show the subhalo-mass cross-correlation considering clumps in different host-halo 
masses. The curves in this case have been normalized to the mean number of clumps in haloes with mass from 10 10 and 10 14 M /h. 



like uU s are the convolution of two NFW profiles, so they 
can also be written analytically (Sheth et al. 2001). 

In Fig. 8, we show the cross-correlation between haloes 
and matter (left panels) as well as between substructures 
and matter (right panels) at the two different redshifts 
z = 0.1 and 2 = 1. These represent the mean redshifts 
of galaxies in the SDSS (Adelman-McCarthy et al. 2006; 
Wang et al. 2007) and in typical weak-lensing surveys, re- 
spectively. In the left panels, we show the halo-matter cross- 
correlation function considering the two components both 
in the 1- and 2-halo terms. We adopt models C3, D2 with 
vine = 0.25 and S2. The solid and dotted curves show the 
total cross-correlation, with distributions SI and S2 for plac- 
ing substructures in the host halo. Notice that, in the host 
mass cross-correlation, the latter choice does not influence 
the total. On the other hand, the cross-correlation between 
substructures and mass is sensitive to the choice of the ra- 
dial subclump distribution. In the right panels, we show the 
cross correlation between clump and mass. Here, too, we 
show the three components of the 1-halo term and the two 
components of the 2-halo term. The solid and dotted curves 
are, respectively, the total predicted signal using models SI 



and S2. From the figure we notice that substructures con- 
tribute a signature on scales of order 30 kpc or so. 

In the left panels of Fig. 9, we show the cross-correlation 
function between haloes and matter, splitting the contri- 
butions into different host-halo mass bins. The integrals in 
Eqs. (38, 39, 40) and (41) have been performed between 
the lower and upper bounds of each bin, while S (k, z) and 
C 1 (k, z) have been integrated over all masses. The solid line 
shows the total power spectrum from haloes in the entire 
mass range, i.e. between 10 10 and lO 14 M /ft. The curves 
are normalized to the mean number of haloes in the entire 
range. 

The long-dashed-dotted curve in the left panels, shows 
the model of Hayashi & White (2008). In their model the 
1-halo term is due to the dark matter density profile of the 
host halo, while the two halo term by the product between 
bias factor and the mass autocorrelation function predicted 
by the linear theory. In their model we consider the mean 
mass predicted by the Sheth & Tormen (1999) mass function 
between 10 10 and 10 14 M Q //i. 

In the right panels, we show the substructure-mass 
cross-correlation considering clumps in the same host-halo 
mass bins. These curves have been normalized by the mean 
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number of clumps in haloes between 10 10 and 10 14 Mq/Ii. In 
this case, showing the predicted contribution from different 
mass bins, we account for halo exclusion in the halo-mass 
cross correlation. An accurate treatment of halo exclusion 
may be numerically tedious. Here, we adopt an approximate 
approach in which the 2-halo term vanishes on scales smaller 
than the virial radius of the mean mass in the bin considered 
(see e.g. (Cacciato et al. 2009) for a similar approach in the 
galaxy-dark matter cross-correlation). Thus, if Mi and M2 
are the lower and upper bounds of the bin, we can write 

1 + £,2 H excl {r\M 1 ,M 2 ) w [1 +bB(r\Mi,Ma)] x W{r\R vir ) 

(48) 

where R V i r is the virial radius of the mean mass between 
Mi and M2 which is used as a typical cut-off scale of the 
2-halo cross-correlation signal here. W(r\R V i r ) is a window 
function equal to 1 if r > R v t r and for r ^ Rvi r - 

The curves in the left panels of Fig. 9 account for this ef- 
fect. The signal near the transition between the 1- and 2-halo 
terms is reduced by 30% due to halo exclusion. The heavy 
dashed line shows the total 2-halo term between 10 10 and 
10 14 Mq/Ji, taking halo exclusion into account. The curve 
tends to zero at a scale corresponding to the virial radius 
of a 10 10 Mq/K halo. The 2-halo term integrals have been 
divided into different bins, for each of which we assign a cut- 
off scale corresponding to the virial radius of a mean halo. 
To estimate the total signal from haloes between 10 10 and 
10 14 M Q /h we sum the contributions from the individual 
mass bins. In the clump-mass cross-correlation, we ignore 
the halo exclusion since the transition between the 1- and 
the 2-halo contributions occurs at larger scales, hence the 
exclusion will not change the signal at all. 



5 SUMMARY AND CONCLUSIONS 

We have shown how the halo-model formalism can be ex- 
tended to account for scatter in halo concentration at fixed 
mass, and for the presence of substructures in dark matter 
haloes. Differences in the mass-concentration relation do af- 
fect the predicted non-linear power spectrum. We quantified 
this by considering three models for the M-c relation, as well 
as a deterministic and a stochastic models for the concen- 
tration at fixed mass. Accounting for substructure means 
that the 1-halo and 2-halo terms should be written as sums 
of four and three types of pairs, respectively. If one uses re- 
alistic models for these different pair-types, then the halo 
model calculation is in reasonable agreement with the non- 
linear power spectrum measured in the Millennium Simula- 
tion, over a range of redshifts. 

The cross-correlation between halos and mass can also 
be estimated using the halo-model formalism. We have 
shown how, also in this case, substructures can be taken 
into account and how the cross-correlation signal between 
clumps and mass can be split into different contributions. 

The agreement with simulations is not yet at the per- 
cent level, so there is room for improvement. The simula- 
tions show an excess at around k ~ 1 relative to our model 
- this is almost certainly not due problems with our sub- 
structure calculation, since it is on larger scales than those 
on which substructures dominate. The discrepancy is also 
unlikely to be due to our decision to ignore 'assembly bias' 
effects. Rather it may be due to the fact that our 2-halo 



term is too crude - on these smaller scales, one must in- 
clude the effects of nonlinear bias, perhaps following Smith 
et al. (2007). 

• Our formulation of the halo model assumes that the 
mass fraction in substructures is a deterministic function of 
the halo concentration: in practice, there is a distribution of 
f s values at fixed c and M. Allowing for this stochasticity 
simply adds one more integral in most of our expressions, 
that can matter at the percent level on the small scales of 
interest here (for essentially the same reason that scatter 
in c and fixed M matters at the ten percent level on small 
scales). 

• We only include some of the effects of the known cor- 
relation between subhalo concentration and distance from 
the center of the host. Specifically, our current implementa- 
tion accounts for the fact that the spatial dependence means 
subhalo concentrations (at fixed mass) are stochastic, but it 
does not include the fact that this is correlated with distance 
from host halo center. Implementing a spatially dependent 
concentration is difficult in Fourier space, so doing so in con- 
figuration is probably the best way forward. 

• We also do not account for the fact that subhalos them- 
selves have substructure. This will matter most for the con- 
tribution which comes from pairs which are in the same sub- 
halo (for the same reason that substructure is not important 
for the 2-halo term). 

• Halo exclusion matters for the cross-correlation signal. 
There are exclusion effects associated with the subhalos too, 
which we currently ignore. 

• Although the sum of the smooth and clumpy compo- 
nents should give an NFW profile, our current implemen- 
tation assumes the smooth component follows and NFW 
profile, whereas the clumpy one does not. Hence there is no 
guarantee that the sum of the two actually is NFW at high 
precision. While this is relatively easy to fix (simply define 
the profile of the smooth component to be the required NFW 
minus the clumpy component), we have not done so. 

• And finally, at this level of precision, halo shapes also 
matter. In particular, the distribution of f s and c at fixed 
M may depend on halo shape - once this has been quanti- 
fied in simulations, it may become necessary to extend our 
formalism to include shapes. 

The formalism presented in this paper can be used in 
subsequent studies to make more realistic predictions for 
the galaxy-galaxy lensing signals and for the lensing conver- 
gence power spectrum. By identifying the different contri- 
butions due to the smooth and clumpy components of the 
dark-matter haloes, this will also allow us to quantify the 
contributions the convergence power spectrum not just by 
different halo masses, but also by different substructures. 
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